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USE OF THE WEST-1 WIND TUR8INE SIMULATOR TO PREDICT 
BLADE FATIGUE LOAD DISTRIBUTION 


David C. Janetzke 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135 


SUMMARY 

The WEST-1 wind turbine simulator Is an Implementation of mathematical 
models for an aeroelastlc horizontal-axis wind turbine rotor using current 
hybrid electronics technology. High-speed digital and analog circuitry enable 
the computation of complex dynamic characteristics and performance of a wind 
turbine In real time. To test the ability of WEST-1 to predict blade fatigue 
load distribution, actual wind signals were fed Into the simulator and the 
response data were recorded and processed In the same manner as actual wind 
turbine data. The WEST-1 simulator was operated In a stable, unattended mode 
for six hours. The probability distribution of the cyclic flatwise bending 
moment for the blade was comparable to that for an actual wind turbine In 
winds with low turbulence. The Input from a stationary anemometer was found 
to be Inadequate for use In the prediction of fatigue load distribution for 
blade design purposes and modifications are necessary. 


INTRODUCTION 

As part of the federal wind energy program, funded by the Department of 
Energy, various analytical tools have been developed for use In the design of 
horizontal-axis wind turbines. One of these developments Is a wind turbine 
simulator called WEST-1, an acronym for Wind Energy System Time Domain-Unit 1 
(ref. 1). The WEST-1 simulator Is a hard-wired hybrid computer specifically 
designed to solve. In real time, the complex linear and nonlinear differential 
equations for the dynamic response of a wind turbine rotor. This real-time 
response permits an analyst to examine the dynamic behavior of a wind turbine 
design with steady state and transient wind conditions for extended periods of 
time. 


The objective of this study Is to evaluate the capability of the WEST-1 
simulator In the prediction of blade fatigue load distribution when measured 
wind data are used for Input. The approach used was to (1) model the DOE/NASA 
Mod-0 rotor on the simulator; (?) use recorded wind speed data from three 
anemometers as Input to the simulator; (3) record, over a six hour period, 
eleven channels of simulator data such as rotor torque and power, blade bend- 
ing moments, wind speed, and wind shear; (4) statistically process the output 
data particularly to determine the probability distribution of the cyclic 
flatwise bending moments at the blade root during the six hour period; and 
(5) compare the blade cyclic flatwise moment probability distribution from the 
simulator with that measured on a comparable wind turbine rotor under field 
test conditions. 
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This paper Includes brief descriptions of the WEST-1 simulator and test 
setup, results from the simulator, and a comparison of blade load distribution 
obtained from the simulator and that obtained from a field-tested wind turbine. 


DESCRIPTION OF WEST-1 

The WEST-1 simulator models an aeroelastlc horizontal-axis wind turbine 
rotor using current hybrid electronic technology (ref. 1). The mathematical 
models are Implemented with a variety of analog, digital, and hybrid computa- 
tional devices. 


Hardware 

The components of the WEST-1 simulator are housed In a printed circuit 
(PC) card cage enclosure shown In figure 1. The enclosure Is 50 cm (20 In.) 
high, 52 cm (21 In.) wide, and 55 cm (22 In.) deep. It holds two card cage 
trays and one power supply tray. The top tray contains circuits which model 
wind speed, wind shear, tower wake and wind retardation. The middle tray con- 
tains all the circuits for the complex aeroelastlc rotor model. The bottom 
tray contains the power supply and a test Instrument for checking and cali- 
brating the various types of PC cards. 

Each PC card contains a number of electronic devices which perform a 
specific mathematical function. Analog devices Include multipliers, summers. 
Integrators, and potentiometers. Digital devices Include logic Inverters, NOR 
gates, J-K fllp-flcps, and one shots. Hyblrd devices Include switches, multi- 
plexors, sample-holds, comparators, and mode-controllable Integrators. 

WEST-1 contains one hundred and thirty-two (132) PC cards and uses over 
sixteen hundred devices. The PC cards are plugged Into a card cage assembly 
called a "pin Diane". One side of the pin plane receives the PC card edge 
connectors and connects each card circuit to a pin on the opposite side. The 
pins are wired to each other with a process called "wire-wrapping" to form the 
desired circuits. 


Mathematical Models 

As presently configured, WEST-1 has the capability to model a horizontal- 
axis wind turbine rotor with the following features: one to four blades; one 

aeroelastlc degree-of-freedom per blade (hinged or cantilevered blade on a 
rigid hub); radial distributions of chord, twist, and mass; full or part-span 
pitch angle; constant rotor speed; linear wind shear; tower wake; and yaw 
angle. An auxiliary analog/digital unit can be used to model the dynamic 
response of other wind turbine components such as the power train, tower, and 
pitch control system. 

The WEST-1 simulator uses a single degree-of-freedom modal representation 
for the aeroelastlc properties of each rotor blade. The first out-of-plane or 
flatwise bending mode of the rotor blade Is used. Aerodynamic and Inertial 
forces are computed with blade element models. These forces are Integrated 
radially to obtain generalized forcing functions and rotor shaft loads. The 
shaft loads from each blade are summed to produce total rotor shaft loads. 
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The air motion models Implemented In the WEST-1 simulator Include trans- 
lational and rotational components of the wind with respect to the rotor, 
tower wake, and retardation of the wind by the rotor. A linear wind shear Is 
modeled by the lateral rotational component of the wind. The tower wake Is 
modeled by an adjustable step change In wind speed over an adjustable sector 
of the rotor disk. The Glauert moment i>:> model Is used to determine the wind 
retardation. 


DESCRIPTION OF TEST SETUP 

The WEST-1 simulator was set up for testing In the Mod-0 Wind Turbine 
Test Facility at the NASA Plum Brook Station In Sandusky, Ohio. Equipment, In 
place at this site, was used to provide wind Input, blade pitch angle control, 
and record the output data. 


Rotor Configuration 

The WEST-1 simulator was programmed, by adjustment of approximately three 
hundred potentiometers, to model the DOE/NASA Mod-0 wind turbine rotor with 
aluminum blades (ref. 2). The basic characteristics of this rotor are: 

(1) Two blades with a rotor diameter of 38 m 

(2) Blade cone angle of 7° 

(3) NACA 230XX series airfoil 

(4) Blade chord distribution of 0.46 to 1.37 m (tlp-to-hub) 

(5) Blade twist angle distribution of -2° to 32* 

(6) Blade thlcxness-to-chord ratio distribution of 0.12 to 0.4*» 

(7) Blade natural frequency of 1.82 Hz for the first flatwise bending mode 

(8) Rotor speed of 40 rpm 

(9) Rated generator power output of 100 kW at wind speed of 8 m/s 


WEST-1 Input Signals 

The wind speed and wind shear Input signals were taken from anemometer 
data at the site which had been previously recorded on analog tape. These 
recorded data provided over six hours of continuous time-varying wind Input. 
Signals from three separate anemometers were used. The overall wind speed 
Input was taken from the signal of an anemometer located at the elevation of 
the Mod-0 wind turbine hub. The difference In the signals from the other two 
anemometers located at elevations of 27 and 59 m provided the linear wind 
shear Input. 

The blade pitch angle Input was controlled by a microprocessor external 
to the WEST-1. The microprocessor was programmed to monitor the power output 
signal and change the pitch angle Input signal as needed to maintain rated 
power output In high wind conditions. The blade pitch angle was held constant 
near zero degrees In low wind conditions. 

Other Input parameters were set by controls on tht front of the simulator. 
The rotor shaft speed was set at a constant 40 rpm. The tower wake was set to 
reduce the wind speed by 35 percent over a 24° sector based on wind tunnel 
tests of the tower model. The yaw angle wa« set at 0°. 
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WEST-1 Output Data Processing 


The data from the simulation were record'd and processed In the same 
manner as those from actual wind turbine tests (ref. 3). Both Input and output 
signals were digitized, screened for maximum and minimum values In each rotor 
revolution, and then recorded on magnetic tape. Recorded Input data Included 
wind speed, wind shear rotor shaft speed, yaw angle, and blade pitch angle. 
Recorded output data Included rotor power, rotor shaft torque, blade root flat- 
wise and chordwlse bending moments, and blade tip deflections. 

Subsequent processing produced midpoint and cyclic parameters of the re- 
corded data according to the following equations: 

midpoint * (maximum ♦ m1n1mum)/2 
cyclic * (maximum - m1n1mum)/2 

The final processing yields a microfiche record containing the time history 
and two statistical summaries of each parameter. The time history plots 
contain the time history of thirty-second averages. One statistical summary 
Is a distribution of the parameter partitioned Into subsets of wind speed. 

The other summary Is a cumulative probability distribution of the parameter. 


RESULTS AND DISCUSSION 

A continuous unattended 6-hour run of the WEST-1 simulator was made with 
pre-recorded anemometer signals used for wind Input and with a microprocessor 
used for blade pitch angle control. The simulator operated In an automatic 
and stable mode throughout the entire period. Data were recorded and pro- 
cessed for analysis and comparison with actual wind turbine results. 


Wind Input Time Histories 

The time histories for the midpoint and cyclic values of the wind speed 
Input to the WEST-1 simulator are shown In figure 2. The midpoint values 
varied continuously and ranged from 5 to 15 m/s. The cyclic values were rela- 
tively constant at a low level of about 0.2 m/s. The time histories for the 
wind shear Input In figure 3 show very similar characteristics. The midpoint 
wind shear values varied frequently between zero and 4 m/s while the cyclic 
values were fairly constant around 0.2 m/s. 

The variability of the midpoint wind speed values Is characteristic of a 
moderately turbulent wind. The low cyclic values are due to normal signal 
noise and are Indicative of a lack of substantial change within one rotor 
revolution at a stationary anemometer. 


Blade Bending Moment Output 

The blade flatwise bending moment output was selected to Indicate the 
capability of WEST-1 In the prediction of blade fatigue loads. It Is re- 
presentative of the primary structural response of a horizontal axis wind 
turbine to variable wind conditions. The aeroelastlc response of the blade 
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first flatwise bending mode Is modeled In WEST-1. The cumulative probability 
distribution of the flatwise bending moments Is of particular Importance In 
the determination of blade fatigue loads. 

Time histories for midpoint and cyclic values of blade flatwise bending 
moments from the WEST-1 run are shown In figure 4. Both curves reflect the 
variability of the wind Input. The partloned distributions of these bending 
moment values with wind speed are given In figure 5. The pair of symbols on 
the vertical bars Indicate the confidence Interval (at the 0.95 level) for the 
differences In the median values for adjacent wind speed subsets. The hori- 
zontal tabs Indicate the 16th and 84th percentile values. 

In figure 5, the midpoint values show a peak near the rated wind speed of 
8 m/s. This reflects the reduction of thrust loads In lower wind speeds at a 
fixed blade pitch angle and In higher winds speeds where the pitch control acts 
to limit the power output. This confirms that the blade pitch angle control 
was functioning properly. 

The cyclic values In figure 5 Increase with wind spac'd and are not Influ- 
enced by the pitch control. This Is reasonable because these cyclic values 
are primarily due to the wind shear and tower wake Input, and the tower wake 
effect Increases with wind speed. 

The cumulative probability distributions for the midpoint and cyclic 
values of the blade flatwise bending moments from the WEST-1 run are given 
In figure 6. The midpoint values range from -77 000 N-m (-57 000 lb-ft) at 
the 0.01 percentile to 24 000 N-m (18 000 lb-ft) at the 99.99 percentile. 

The cyclic values range from 25 000 N-m (19 000 lb-ft) to 78 000 N-m 
(58 000 lb-ft). 


Comparison of Blade Load Distributions 

The wind turbine test data, used here for comparison with the WEST-1 
simulator results, were measured on the DOE/NASA Mod-OA wind turbine at 
Clayton, New Mexico. An extensive set of data was recorded and processed and 
has been reported In reference 4. Although rated at 200 kW, the Nod-OA rotor 
Is essentially the same as the Mod-0 rotcr configuration modeled In the WEST-1 
simulator. As noted In the discussion of figure 5, this difference In power 
control does not slgnflcantly affect blade flatwise cyclic loads. 

A normalized plot of WEST-1 data and two sets of Mod-OA test data are 
presented In figure 7. The probability distributions for the blade cyclic 
flatwise bending moments of these sets are normalized with respect to their 
median (50th percentile) values. This normalizing removes differences due to 
mean wind speed, wind shear, pitch angle, and radial location, permitting 
direct comparison of the distribution slopes. The slopes can be used to cal- 
culate the log standard deviations for the various probability distributions. 

The WEST-1 data set, shown In figure 7, Is the normalized version of that 
presented i figure 6. The median wind speed was 10.0 m/s for the six-hour 
test run. The turbulence level was moderate. 
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The first Mod-OA data set In figure 7 Is a comprehensive set compiled 
from 43 hours of test data (ref. 4). It Includes data from tests covering the 
full range of wind speeds and turbulence levels. The median wind speed was 
8.0 m/s for this data set. All levels of turbulence (low, moderate, and high) 
were encountered In these tests. 

The second Mod-OA data set In figure 7 Is a subset of the first. These 
data were obtained In a six-hour test run during which the wind turbulence 
level was low. The median wind speed was 7.4 m/s. 

Straight-line curve fits through the data In figure 7 Indicate log-normal 
distributions for the blade cyclic flatwise bend'ng moments. The slope of the 
WEST-1 curve fit Indicates a log standard deviation of 0.077. The comprehen- 
sive Mod-OA dataset has a significantly higher log standard deviation of 0.201 
The Mod-OA subset with low-turbulence-level wind conditions has a log standard 
deviation of 0.073 comparable to that of the WEST-1. 

The blade fatigue load distribution produced by the WEST-1 simulator Is 
unconservative with regard to both Mod-OA data sets. The WEST-1 distribution 
Is comparable to the six-hour Mod-OA distribution despite the lower wind speed 
and low turbulence level of the Mod-OA run. The blade load distribution of 
the comprehensive Mod-OA data set Is substantially wider than that of the 
WEST-1 because of the wide range of wind speed and turbulence levels. However 
It Is this wider range of loads which the simulator must predict In order to 
become a useful design tool. 

Apparently, based on the wind Input and blade load response, the wind 
Input derived from three stationary anemometers does not adequately model the 
character of the wind across the rotor disk. The low cyclic values of the 
wind Input from a moderately turbulent wind produce cyclic blade loads com- 
parable only to actual loads In low turbulence conditions. A wind character- 
istics study has shown that turbulence encountered by a point on a rotating 
blade Is quite different than turbulence measured by a stationary anemometer 
(ref. 5). In that study, the spectral analysis of rotatlonally sampled wind 
data Indicates significant amounts of energy In turbulence fluctuations at 
multiples of the rotational speed. These characteristics are not found In 
stationary anemometer data and are evidently needed In the wind Input for the 
simulation of blade response to turbulent wind conditions. Thus modifications 
to the wind Input are necessary to produce a blade fatigue load distribution 
useful for design purposes. 


CONCLUSIONS AND RECOMMENDATIONS 

The following conclusions are made from the results of the WEST-1 simu- 
lator run and the comparison with actual wind turbine data: 

(1) The WEST-i wind turbine simulator can operate In a stable unattended 
mode continuously for extended periods of time with time-varying wind Input 
and microprocessor blade pitch angle control. 

(2) The WEST-1 simulator with stationary anemometer wind speed Input 
produces a blade fatigue load spectrum similar to wind turbine test data for 
low wind turbulence conditions. 


6 



(3) The wind signals from stationary anemometers are not adequate for 
Input to the WEST-1 In producing blade fatigue load distribution for use In 
blade design. 

It Is recommended that further testing of the WEST-1 simulator be made 
with modifications to the wind speed Input which would expand the range and 
Increase the variations within each rotor revolution. 
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Figure 3. - Wind shear input to WEST-1 simulator 
from difference of anemometer signals at eleva- 
tions of 59 and 27 meters. 



Figure 4 - Time history of blade flatwise bending moments 
from WEST-1 simulator. 
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Figure 5. - Partitioned distribution of blade flatwise bending 
moments with wind speed from WEST-1 simulator. 
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Figure 6. - Cumulative probability distribution of blade flat- 
wise bending moments from WEST-1 simulator. 
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Figure 7. - Normalized blade cyclic flatwise bending 
moment probability distribution tor WEST-1 simula- 
tor and comparable Mod-OA wind turbine test data. 



